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We suggest a novel conceptual approach for describing the properties of waveguides and circuits in 
photonic crystals, based on the effective discrete equations that include the long-range interaction 
effects. We demonstrate, on the example of sharp waveguide bends, that our approach is very 
effective and accurate for the study of bound states and transmission spectra of the photonic-crystal 
circuits, and disclose the importance of evanescent modes in their properties. 



One of the most promising applications of photonic 
crystals is a possibility to create compact integrated op- 
tical devices [jj] , which would be analogous to integrated 
circuits in electronics, but operating entirely with light. 

Usually, the properties of photonic crystals and 
photonic-crystal waveguides are studied by solving 
Maxwell's equations numerically, and such calculations 
are time consuming. Moreover, the numerical approach 
does not always provide a good physical insight. The 
purpose of this Letter is to suggest a novel approach, 
based on the effective discrete equations, for describing 
many of the properties of the photonic-crystal waveguides 
and circuits, including the example of the transmission 
spectra of sharp waveguide bends. The effective discrete 
equations we derive below are somewhat analogous to the 
Kirchhoff equations for electric circuits. However, in con- 
trast to electronics, in photonic crystals both diffraction 
and interference become important, and thus the result- 
ing equations involve the long-range interaction effects. 

To introduce our approach, we consider a two- 
dimensional (2D) photonic crystal consisting of infinitely 
long dielectric rods arranged in the form of a square 
lattice with the lattice spacing a. We study the light 
propagation in the plane normal to the rods, assuming 
that the rods have a radius ro = 0.18a and the dielec- 
tric constant Sq — 11.56 (this corresponds to GaAs or 
Si at the wavelength ~ 1.55 firn). For the electric field 



E(x,t) 



E(x\lo) polarized parallel to the rods, 



Maxwell's equations reduce to the eigenvalue problem 



E(x\to) = 



(1) 



which can be solved by the plane- wave method (2). A 
perfect photonic crystal of this type possesses a large 
(38%) complete band gap (between uj — 0.303 x 2irc/a 
and uj = 0.444 x 2irc/a), and it has been extensively 
employed during last few years for the study of bound 
states in waveguides and bends |J , transmission of light 
through sharp bends , branches Q and channel drop 
filters juj , nonlinear localized modes in straight waveg- 
uides || and perfect photonic crystals ||. Recently, this 
type of photonic crystal with a 90° bent waveguide was 
fabricated experimentally in macro-porous silicon with 



a = 0.57 fim and a complete band gap at 1.55 fim fl(i[| . 

To create a waveguide circuit, we introduce a system 
of defects and assume, for simplicity, that the defects are 
identical rods of the radius r d (with Ed) located at the 
points x m , where m is the index number of the defect 
rods. In the photonic crystal with defects the dielectric 
constant e(x) can be presented as a sum of the periodic 
and defect-induced terms, i.e. s(x) = e p (x) + Sd(x) , and, 
therefore, Eq. (|l]) can be written in an integral form 

E{x\u) = f d 2 y G{x,y\uj)e d {y)E{y\u) , (2) 



where G(x,y\ui) is the Green function (see, e.g., S). 

The integral equation (|^) can be solved numerically 
in the case of a small number of the defect rods. How- 
ever, such calculations become severely restricted by the 
current computer facilities as soon as we increase the 
number of the defect rods in order to create photonic- 
crystal waveguides, waveguide bends, and branches [Q-Q. 
Therefore, our primary goal in this Letter is to develop 
a new approximate physical model that would allow the 
application of fast numerical techniques combined with a 
reasonable accuracy and the further possibility to study 
nonlinear photonic crystals and waveguides. 

When the defects support monopole modes, a reason- 
ably accurate model can be derived by assuming that 
the electric field inside a defect rod remains constant. In 
this case, we can average the electric field in the inte- 
gral equation (0) over the cross-section of the rods ||[ll) , 
and derive an approximate matrix equation for the am- 
plitudes of the electric field E n (uj) = E(x n \cu) at the 
defect sites, 

Y,M nim (w)E m = , 

m 

M„ im (w) = e d J„,m(w) - 8 n>m , (3) 
where <5„. m is the Dirac's delta function, and 

J n ,m(v) = (j^j j d 2 y G(x n ,x m + y\uj) (4) 



1 



is a coupling constant determined through the Green 
function of a perfect 2D photonic crystal MM . 
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FIG. 1. 
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Dispersion relation for the 2D photonic-crystal 



waveguide (shown in the inset) calculated by the super-cell 
method [g| (dashed), and from the approximate equations 
(§)-(§) for L=7 (solid) and L=l (dotted). The gray areas 
are the projected band structure of a perfect 2D crystal. 



To check the accuracy of the approximate model (^), 
first we consider a single defect located at the point xq. 
In this case, Eq. (J|) yields Jo,o(^d) = l/ e <2' an d this 
expression defines the frequency iOd of the defect mode. 
For example, applying this approach to the case when we 
have a defect created by a single removed rod, we obtain 
the frequency uJd = 0.391 x 2irc/a which differs by only 
1% from the value u>d — 0.387 x 2nc/a calculated with 
the help of the MIT Photonic-Bands numerical code Q . 

A single-mode waveguide can be created by removing a 
row of rods (see the inset in Fig. [l]). Assuming that the 
waveguide is straight (M n , m = M„_ m ) and neglecting 
the coupling between asunder defect rods (i.e. M„_ m = 
for all | n — to | > L), we rewrite Eq. (^) in the transfer- 
matrix form, F n+ i = TF n , where we introduce the vector 
F n = { E n , E n —i, ... , £"ri-2L+i } and the transfer matrix 
T = {Ti j} with the non-zero elements 

Tx.-H = -%# for j = l,2,...,2L, 



T 



M l {lo) 
1 for j = 1,2,...,2L- 1 



Solving the eigenvalue problem 

f{oj)& = exp{ik p (uj)} & , 



(5) 



(6) 



we can find the 2L eigenmodes of the photonic-crystal 
waveguide. The eigenmodes with real wavenumbers 
k p (u>) correspond to the propagating waveguide modes. 
In the waveguide shown in Fig. Jl], there exist only two 
such modes (we denote them as < I> 1 and $ 2 ), propagating 
in the opposite directions (fei = — fc 2 > 0). In Fig. [I] we 
plot the dispersion relation ki(u>) found from Eq. (g) for 
the nearest-neighbor interaction (L=l) and also taking 
into account interaction between seven neighbors (X=7); 
we compare the results with those calculated directly by 



the super-cell method [||. As soon as we go beyond the 
approximation of the nearest neighbors and take into 
account the coupling between several defect rods, Eqs. 
(H) - © provide very accurate results for the dispersion 
characteristics of the photonic-crystal waveguides. We 
verify that this conclusion is also valid for multi-mode 
waveguides, e.g. those created by removing several rows 
of rods. 
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FIG. 2. Electric field E n for two bound states supported 
by a 90° waveguide bend (shown in the top). Center of the 
bend is located at n = 0. 



In addition to the propagating guided modes, in 
photonic-crystal waveguides there always exist evanes- 
cent modes with imaginary k p . These modes, which can- 
not be accounted for in the framework of the nearest- 
neighbor approximation, remain somewhat "hidden" in 
straight waveguides, but they become important in 
more elaborated structures such as waveguide bends and 
branches. Importantly, our model does take into account 
all such effects. 

We consider the simplest case of a waveguide bend, 
where the evanescent modes manifest themselves in two 
different ways. First of all, they create localized bound 
states in the vicinity of the bend. As was shown in Ref. 
H , in the cases when the waveguide bend can be consid- 
ered as a finite section of a waveguide of different type, 
the bound states correspond closely to cavity modes ex- 
cited in this finite section. However, such a simplified 
one-dimensional model does not describe correctly more 
complicated cases, even the bent waveguide depicted in 
Fig. | H . The situation becomes even more complicated 
for the waveguide branches In contrast, solving Eq. 
(Q) we can find the frequencies and profiles of the bound 
states excited in an arbitrary complex set of defects. As 
an example, in Fig. ^| we plot the profiles of two bound 
states (cf. Fig. 9 in Ref. |H|). The frequencies of the 
modes are found from Eq. (||) with the accuracy of 1.5%. 

Additionally, the evanescent modes determine the non- 
trivial transmission properties of the waveguide bends 
which can also be calculated with the help of our dis- 
crete equations. To demonstrate this, we consider a bent 
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waveguide consisting of two coupled semi-infinite straight 
waveguides with a finite section (an arbitrary complex set 
of defects) between them. The finite section includes a 
bend with a safety margin of the straight waveguide at 
both ends. We assume that the defect rods inside this 
segment are characterized by the index that runs from a 
to b, and the amplitudes E m (m = a, ...,&) of the elec- 
tric field near the sites of the removed rods are all un- 
known. We number the guided modes (^J) in the following 
way: p = 1 corresponds to the mode propagating in the 
direction of the waveguide bend (for both ends of the 
waveguide), p — 2 corresponds to the mode, propagating 
in the opposite direction, p = 3, L + 1 correspond to 
the evanescent modes which grow in the direction of the 
bend, and p = L + 2, 2L correspond to the evanescent 
modes which decay in the direction of the bend. Then, 
we can write the incoming and outcoming waves in the 
semi-infinite waveguide sections as a superposition of the 
guided modes: 

L + l 

E in ,£,1 , $ 2 , A m $ p f 7 x 

^ a— in > ' ^ a — m ~ / t yy p ^ a — m > V 1 / 

for m = a — 2L, a — 1, and 

p=3 

for m = b+ 1, b + 2L, where A™ and \° ut are unknown 
amplitudes of the evanescent modes growing in the di- 
rection of the bend, and t and r are unknown amplitudes 
of the transmitted and reflected propagating waves. We 
take into account that the evanescent modes growing in 
the direction from the bend vanish, and assume that the 
amplitude of the incoming plane wave $ 1 is normalized 
to the unity. Now, substituting Eqs. (@)-(§) into Eq. 
we obtain a system of linear equations with 2L + b — a + 1 
unknown. Solving this system, we find the transmission 
\t\ 2 and reflection |r| 2 coefficients. 

In Fig. [j] we present our results for the transmission 
spectra of several types of bent waveguides, as in Ref. [Q, 
where the possibility of high transmission through sharp 
bends in photonic-crystal waveguides was first demon- 
strated. As is clearly seen, Eqs. (^)-(^) provide a very 
accurate method for calculating the transmission spectra 
of the waveguide bends. 

In conclusion, we have suggested a novel conceptual ap- 
proach for describing the properties of photonic-crystal 
waveguides and circuits, including the transmission spec- 
tra of sharp bends. The effective discrete equations we 
have introduced here emphasize the important role of the 
evanescent modes in the photonic-crystal circuits, and 
they can be applied to study more complicated prob- 
lems such as transmission in waveguide branches, chan- 
nel drop filters, nonlinear localized modes in nonlinear 
waveguides, and so on. 

The authors are indebted to S.H. Fan for useful com- 
ments and to A. Mekis for providing the data from 



Ref. pi]. The work has been partially supported by the 
Australian Research Council. 
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FIG. 3. Reflection coefficients calculated by the fi- 
nite-difference time-domain method (dashed, from Ref. El) 
and from Eqs. (||)-(||) with L = 7 (full lines) and L = 1 (dot- 
ted, only in the top plot), for four different bend geometries. 
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